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Abstract 

The  results  of  an  experiment  examining  the  forward  propagation  of  an  acoustic 
signal  through  a  buoyant  plume  are  discussed.  Two  distinct  testing  sights  were  used.  One 
made  use  of  a  small  fresh  water  tank  in  NUWC  to  provide  a  controlled  plume.  The  other 
used  a  larger  salt  water  tank  at  WHOI  to  create  a  more  realistic  oceanic  model.  Using  the 
Bom  and  Rytov  approximations,  an  estimation  of  the  effects  of  the  laminar  plume  on  the 
propagated  signal  are  shown.  As  the  plume  moves  from  laminar  to  turbulent,  the 
scintillation  index  and  the  Fourier  transform  of  the  magnitude  square  response  provide 
insight  into  the  nature  of  the  transition.  Finally,  from  the  turbulent  response  a  model  for 
the  scattering  function  is  developed. 
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Chapter  1 


INTRODUCTION 

1.1  Statement  of  Problem 

This  thesis  examines  the  effect  on  a  pulsed  continuous  wave  (CW)  signal  of 
thermally  generated  laminar  and  turbulent  plumes.  Initially  the  forward  scattering  through 
a  laminar  plume  is  modeled  using  the  Bom  and  Rytov  approximation.  Then  a  model  using 
scattering  from  a  homogeneous  cylinder  is  used  to  allow  easy  comparisons  for  different 
thermal  profiles  as  the  plume  moves  from  laminar  to  turbulent.  The  Fourier  transform  of 
the  time  sequence  of  the  magnitude  response  squared  of  the  signal  train  is  used  to 
investigate  the  effective  amplitude  modulation  of  the  index  of  refraction.  Then,  the 
scintillation  index  is  used  to  discuss  the  point  at  which  the  plume  becomes  turbulent. 
Finally,  the  scintillation  index  and  Fourier  transform  of  the  time  sequence  magnitude  are 
used  to  compare  data  obtained  under  highly  divergent  conditions. 

1.2  Overview  of  Acoustic  Turbulent  Studies 

The  genesis  of  this  research  was  a  discussion  concerning  a  low  frequency 
modulation  of  a  shallow  water  telemetry  signals.  The  problem  was  resolved,  but  the 
question  of  the  effects  of  buoyant  plumes  on  line  of  sight  propagation  had  been  raised. 
Ocean  buoyant  plumes  can  arise  from  a  freshwater  spring,  a  heat  vent,  or  even  the  wake  of 
body  moving  through  a  layered  medium.  These  sources  vary  considerably  depending  on 
their  extent;  the  degree  to  which  they  alter  the  effective  index  of  refraction;  and  the  amount 
of  entrained  particulate  scatterers.  In  recent  years  a  significant  amount  of  research  has 
concerned  observations  of  hydrothermal  vent  fields.  However,  available  literature  does  not 
indicate  as  intensive  a  research  into  other  buoyant  plumes.  The  investigation  of  these 
phenomena  may  be  enhanced  from  the  results  obuuned  in  this  thesis. 
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Aside  from  the  observation  of  the  hydrothermal  iields,  most  research  on  turbulent 


fields  has  focused  on  propagation  through  turbulent  media.  Extensive  study  of 
elecU’omagnetic  wave  propagation  through  turbulent  atmosphere  exists. *  Chernov-  and 
Tatarski^  provide  the  basics  of  wave  interaction  with  turbulent  and  random  media.  Each 
takes  a  different  approach  to  similar  subjects  and  if  one  proves  difficult  to  understand  the 
other  will  provide  the  necessary  insight  Acoustic  propagation  has  been  treated  by  Wenzel 
and  Keller^  and  more  recendy  by  Yang  and  McDaniel,^  but,  in  general  the  electromagnetic 
literature  has  been  sited  and  used  where  appropriate.  Flattd  discussed  the  appropriate 
regions  for  these  applications^  to  ocean  acoustics. 

The  work  of  Mintzer  in  1953  provides  a  notable  exception  to  the  application  of 
electromagnetics  direedy  to  ocean  phenomena.  Mintzer”^  attempted  to  provide  a  model  for 
the  results  of  experimental  work  of  Sheehy  who  had  used  a  series  of  sound  pulses  to 
investigate  the  variaUon  in  the  propagadon  of  sound  in  a  inhomogeneous  media  for  the 
Office  of  Naval  Research  in  1948.  The  work  is  notable  not  just  for  the  use  of  pulsed  waves 
but  also  for  the  use  of  "the  coefficient  of  variadon;"  which  has  the  exact  form  of  the 
scindllation  index  except  for  the  use  of  the  absolute  value  of  the  pressure  field  ( |pj ) 

instead  of  the  intensity  ( I ). 


^  Akira  Ishimaru,  Wave  Propagation  and  Scattering  In  Random  Media,  Vol .  2 
(New  York, San  Francisco,  London:  Academic  Press, 1978) , Part  IV. 

^Lev  A.  Chernov,  Wave  Propagation  in  A  Random  Medium,  (New  York, 

Toronto,  London:  McGraw-Hill,  1960) . 

^V.I.  Tatarski,  Wave  Propagation  in  a  Turbulent  Medium,  (New  York, 
Toronto,  London:  McGraw-Hill,  1961) . 

Alan  R.  Wenzel  and  Joseph  B.  Keller,  "Propagation  of  Acoustic  Waves  in 
a  Turbulent  Medium, "  The  Journal  of  the  Acoustical  Society  of 
America, 50, No.  3  (part  2) (January  1971) , 911-920  . 

^  J.H.  Tarng,  C.C.  Yang,  and  S.T.  McDaniel , "Acoustic  Beam  Propagation  in 
a  Turbulent  Ocean, "The  Journal  of  the  Acoustical  Society  of  America,  84, 
No.  5  (November  1988) , 1808-1812 . 

®  Stanley  M.  Flatte, "Wave  Propagation  Through  Random  Media: 

Contributions  From  Ocean  Acoustics,"  Proceedings  of  The  IEEE,  71, No.  11 
(November  1983)  1267-1294. 

^  David  Mintzer,  "Wave  Propagation  in  a  Randomly  Inhomogeneous 
Medium . I ," The  Journal  of  the  Acoustical  Society  of  America, 25, No. 5 
(September  1953)  922-927. 
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(1) 


^  (ipr)-i(p) 

i(5>r 

Unfortunately,  only  propagation  through  the  entire  media  was  considered  with  no 
accounting  for  discrete  phenomena.  Forward  scattering  in  the  ocean  through  weak 
fluctuations  is  also  addressed  by  Duda^  with  heavy  emphasis  on  the  scintillation  index. 
Both  Duda  and  Mintzer  look  at  the  increase  in  the  scintillation  index  (coefficient  of 
variation)  with  range  through  the  turbulent  medium. 

There  have  been  several  papers  that  deal  directly  with  turbulence  or  ocean 
turbulence  as  a  local  phenomena.  Rasmussen^  discusses  the  localization  of  turbulent 
phenomena  based  on  the  reverberation  spectrum  spreading  due  to  turbulent  scatterer  motion 
in  contiguous  reverberation  volumes.  The  spreading  of  the  spectra  discussed  involves 
motion  of  the  scatterers  on  the  order  of  tens  of  cm/s  in  order  to  create  a  sufficient  Doppler 
frequency.  Artel'nyi  and  Raevskii'®  look  at  the  effect  of  multiple  scattering  of  ocean 
turbulence  by  looking  at  the  effect  of  small  fluctuations  in  the  sound-velocity-profiJe. 

They  ignore  the  temporal  fluctuations  of  the  turbulence  "since  the  wave  frequency  O)  is 
usually  much  greater  than  the  characteristic  frequencies  of  the  turbulent  fluctuations."** 
Instead  Artel'nyi  and  Raevskii  concentrate  on  the  effects  in  the  modal  spectra. 

Papers  directly  addressing  the  effect  of  plumes  and  turbulent  plumes  on  acoustic 
propagation  have  been  becoming  increasingly  frequent  in  recent  years.  Among  the  earliest 
examples  of  the  work  in  this  field  was  "The  Scattering  of  Sound  in  a  Turbulent  Medium" 


®  Timothy  F.  Duda, "Modeling  Weak  Fluctuations  of  Undersea  Telemetry 
Signals,"  IEEE  Journal  of  Oceanic  Engineering,  16,  No.  1  (January 
1991),  3-11. 

^  Robert  A.  Rasmussen, "  Remote  Detection  of  Turbulence  From  Observations 
of  Reverberation  Spectra, "The  Journal  of  the  Acoustical  Society  of 
America, 63 , Ho .  1  (January  1978)  101-110. 

V.V.  Artel'nyi  and  M.A.  Raevskii, "Multiple  Scattering  of  Low- 
frequency  Sound  Waves  by  Ocean  Turbulence,"  Soviet  Phys .  Acoust.,33,  No. 
1  ( January-February  1987),  4-7. 

Artel’nyi,  p.4. 
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by  Kraichnan*2  which  used  turbulent  velocity  field  confined  to  a  finite  region  and  modeled 
the  effect  using  an  isotropic  field  Ihe  resulting  expression  for  the  scattering  by  isou-opic 
turbulence  is: 


(P(n,  e  ))  = 


TtO) 


sin^  20 


4  sin' 


f-1 

\2) 


+  5' 


to' ^4  sin'  ^  +  5'  j 


,  0)  ,5 


(2) 


and,  demonstrates  that  for  isotropic  turbulence  the  mean  scattered  field  in  the  forward 

direction  is  zero  (0).  F  is  the  ensemble  average  of  the  scalar  of  the  frequency  wavenumber 
response  of  the  field:  F  =  ^|f(k,a))|^^  and  n  is  the  unit  vector  in  a  direction  from  the 

scatterer.  Buoyant  plumes  may  be  internally  isotropic,  but  by  their  very  definition  as 
buoyant,  imply  a  variation  between  the  mean  internal  buoyancy  and  the  surrounding  fluid 
buoyancy.  The  mean  forward  scattered  signal  should  contain  at  least  a  small  variation  in 
the  mean  field  in  the  forward  direction  for  a  buoyant  plume. 

The  off  axis  scattering  of  a  turbulent  plume  for  forced  flow  systems  has  been 
investigated  for  flow  rates  near  650  cm/s  by  Korman*^  and  Beyer.  Korman  used  a 
modified  form  of  (1)  to  estimate  the  received  spectrums.  The  back  scattered  signal  from  a 
particle  laden  turbulent  plume  has  been  investigated  by  Hay*'*.  Hay  arrived  at  a  Gau.ssian 
particle  distribution  for  the  entrained  particulates.  The  Gaussian  distribution  appears  to  be  a 
dominating  form  for  many  turbulent  parameters.  Hay’s  work  however  concentrated  on  the 
backscattered  signal  from  varying  particle  types  and  concentration. 


Robert  H.  Kraichnan, "The  Scattering  of  Sound  in  a  Turbulent  Medium, " 
The  Journal  of  the  Acoustical  "ociety  of  America,  25,  No.  6  (November 
1953),  1096-1104. 

M.S.  Korman  and  R.  T.  Beyer, "The  Scattering  of  Sound  by  Turbulence  in 
Water,  "The  Journal  of  the  Acoustical  Society  of  America, 61,  No.  6  (June 
1980),  1980-1987. 

Alex  E.  Hay, ‘Sound  Scattering  from  a  Particle  Laden,  Turbulent 
Jet, "The  Journal  of  the  Acoustical  Society  of  America, 90,  No.  4  (October 
1991),  2055-2073. 
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Wittek‘5  looked  at  the  amplitude  modulation  created  by  turbulent  free  flow  jet  at 
various  velocities  in  the  forward  scattered  direction.  By  looking  at  the  frequency  spectra  of 
the  forward  scattered  field,  an  estimate  of  the  nature  of  the  turbulence  could  be  obtained. 
The  jet  flow  velocities  used  were  at  a  minimum  of  2.3  m/s,  almost  one  hundred  times  the 
velocity  obtained  in  this  buoyant  plume  experiment  The  frequency  spectrum  obtained  for 
the  forward  scattered  field  intensity  closely  parallels  the  results  expected  in  Ishimaru*^  for  a 
Gaussian  velocity  distribution  with  a  Doppler  shift  due  to  the  mean  velocity  and  frequency 
spreading  due  to  velocity  fluctuation. 

Goodman  and  Oeschger*^  investigated  scattering  from  a  buoyant  laminar  plume  in 
the  off  axis  directions.  Further  investigation  by  Oeschger**  followed.  The  temperature 
profile  of  a  laminar  plume  was  measured  and  verified  to  be  Gaussian.  Then,  using  the  Bom 
approximation,  the  scattered  field  was  determined  in  terms  of  the  Bragg  wave  number  and 
verified  by  experimental  data.  The  researchers  did  not  investigate  forward  scattering  of  the 
signal  through  the  laminar  plume.  In  addition  they  did  not  look  at  the  plume  when  turbulent 
or  in  transition.  The  laboratory  and  equipment  used  in  the  Goodman  and  Oeschger 
experiments  were  also  used  in  the  experiment  for  this  thesis. 

An  understanding  of  laminar  and  turbulent  plumes  on  acoustic  signals  can  be  seen, 
by  using  the  results  for  the  laminar  plume  obtained  by  Goodman  and  Oeschger  and 
extending  them  to  varying  power  levels.  Ultimately,  a  further  understanding  of  the  nature 
of  the  transition  from  laminar  to  turbulent  in  terms  of  their  acoustic  effects  can  also  be 
obtained. 


G.  Wittek, "Amplitudenmodulation  Von  Ultraschallwellen  durch 
turbulence  Wasserstromungen"  Acustica,  27  (1972),  7-14. 

Ishimaru,  p.341. 

L.  Goodman,  J.  Oeschger,  and  D.  Szargowicz, "Ocean  Acoustics 
Turbulence  Study:  Acoustic  Scattering  from  a  Buoyant  Axisymnietric 
Plume, "The  Journal  of  the  Acoustical  Society  of  Ameri ca , 9 1 , No .  6  (June 
1992),  3212-3227. 

John  W.  Oeschger , Acoustic  Scattering  From  a  Laminar  Plume, University 
of  Rhode  Island, 1993 . 
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Chapter  2 


Experimental  Description 

The  effect  of  turbulence  on  the  propagation  of  sound  is  primarily  due  to  the 
variation  in  the  effective  index  of  refraction  across  the  plume.  For  this  experiment  a  heating 
element  is  used  to  create  an  initially  laminar  plume  which  becomes  turbulent  as  power  to  the 
plume  is  increased.  An  acoustic  signal  is  sent  across  the  plume  to  study  the  effects  of  the 
plume  in  the  forward  scattering  direction. 

The  primary  mechanism  for  the  variation  in  the  index  of  refraction  is  the  variation  in 
the  temperature  across  the  plume.  From  measurements  conducted  at  60 v  (1 1.7  Watts),  the 
laminar  plume  was  found  to  have  the  expected  Gaussian  temperature  distribution  within  a  2 
cm  diameter.!^  Once  the  plume  is  turbulent  the  expected  average  thermal  profile  is  also 
Gaussian.20  Using  these  assumptions  the  variation  in  the  laminar  plume  with  power  can  be 
modeled  as  a  cylinder  with  an  effective  internal  index  of  refraction. 

The  turbulent  plume  creates  an  effective  amplitude  modulation  of  the  received  signal 
about  a  mean  value.  The  extent  and  form  of  the  turbulence  can  be  determined  from  the 
frequency  spectrum  of  the  variations  in  the  received  signal  amplitudes.  By  using  the 
scintillation  index  and  the  frequency  spectrum  as  measures  of  the  variations  within  the 
turbulence,  a  rough  estimate  of  the  point  at  which  the  plume  becomes  as  turbulent  as  is 
possible  can  be  determined. 

Two  different  laboratory  set  ups  were  used  to  investigate  the  phenomena.  One 
apparatus  used  the  small  fresh  water  tank  used  by  Oeschger  and  Goodman;  this  provided  a 
controlled  result  with  a  known  thermal  plume.  The  other  used  a  large  1 1.65  m^  salt  water 

Oeschger,  p.  36. 

H.  Rouse, C.  Yih,  and  H.W.  Humphreys, "Gravitaional  Convections  from  a 
Boundry  Source ."Tel lus, A  (1952),  201-210. 
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tank  which  was  allowed  to  stratify  and  develop  internal  waves.  The  internal  waves  were 
generated  by  location  in  the  room  and  internal  and  external  weather  conditions.  The  latter 
tank  was  used  to  investigate  the  robustness  of  the  measures  of  the  turbulence. 

The  two  experimental  systems  completed  essentially  the  same  experiment.  The  face 
to  face  distance  between  the  transducers  is  the  only  experimental  parameter,  outside  of  the 
water  in  the  tank,  that  is  varied.  Both  experimental  systems  allowed  the  transducers  to  be 
placed  face  to  face  at  specific  heights  above  the  heater  element  Because  the  phenomena  is 
discrete  and  in  the  far  field  of  the  transducers  the  experiments  can  be  compared  directly. 

2.1  Fresh  Water  System 

A  tank  made  of  2.5  cm  Plexiglas  63.5  cm  wide  1 18. 1 1  cm  long  and  8 1 .28  cm  deep 
of  nominal  depth  68  cm  with  an  open  top  to  allow  access  to  equipment  was  filled  with 
filtered  tap  water.  Initially  filtering  used  a  Bamstead  nanopure  water  filtering  system  in  the 
gravity  feed  mode  producing  a  1.5  liter  per  minute  flow  rate  while  eliminating  panicles 
down  to  0.2  microns.  Eventually,  a  new  filtration  system  by  Pureflow  with  a  1  micron 
filtration  at  150  liters  per  minute  replaced  the  old  system.  The  new  filtration  system  was  in 
place  for  the  data  set  at  a  height  of  15  cm.  The  systems  were  turned  off  during  the 
experiments  to  prevent  distortions  in  the  plume  due  to  filU'ation  flow.  The  continual  flow 
helped  limit  stratification  of  the  tank  as  well  as  removing  particulates. 

An  anodized  aluminum  frame  with  an  open  internal  ring,  allowing  face  to  face 
transducer  ranges  of  50.96  cm,  was  placed  in  the  tank.  Threaded  stainless  steel  rods 
attached  the  ring  to  the  base.  The  rods  allowed  variation  in  the  transducer  height  above  the 
plume.  The  base  held  the  heating  element  in  the  center  of  the  ring.  The  frame  could  be 
removed  from  the  tank  to  make  adjustments.  The  heating  element,  0.95  cm  diameter  by 
2.5  cm  in  height,  was  rated  at  50  watts  at  120  volts  rms.  The  analog  variable  power  supply 
allowed  variation  from  0  to  140  volts.  By  varying  the  power  to  the  heater  the  plume  varies. 
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Figure  2.2 

The  nature  of  the  plume  ( laminar,  unstable,  or  turbulent )  needed  to  be  determined 
easily  in  order  to  verify  the  acoustic  results.  A  laser  shadow  graph  system  provided  real 
time  qualitative  information  on  the  nature  of  the  plume.  A  5  milliwatt  diode  laser  having  a 
wavelength  of  670  nm  provided  the  light  source.  The  laser  beam  was  expanded  and 
directed  through  the  plume  perpendicular  to  the  acoustic  axis  and  the  plume  axis.  On  the 
opposite  side  of  the  tank  the  beam  was  directed  to  a  screen  from  which  a  video  image  was 
recorded,  digitized,  and  displayed  using  a  Raster  Ops  364  Color  Video  board.  The 
digitized  signal  was  saved  periodically  for  latter  analysis. 

The  transducers  used  in  the  experiment  were  3/4"  Panametrics  transducers  with  1 
Mhz  center  frequencies  with  a  -10  dB  bandwidth  of  400  kHz  to  1.4  MHz.  "Beam  pattern 
data  taken  on  the  transducers  showed  that  they  could  be  well  modeled  by  baffled 
pistons."^*  Using  a  model  of  the  transducer  as  a  plane  circular  piston  the  "U.S.  Institute 
(now  American  National  Standards  Institute)  arbitrarily  decided  that  the  criterion  for  a 


Oeschger,  p.29. 
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circular  piston  will  be  that  the  far  field  begins  at  the  critical  range  Rc  =pa2/l  "22  for 
a=transducer  radius.  This  corresponds  to  a  value  of  19  cm.  Another  approach  for 
determining  the  point  at  which  the  far  field  begins,  uses  the  known  expression  for  the  on 
axis  extrema:  rni=(a2/ml)-(ml/4),  and  then  defines  the  far-field  region  as  the  region  beyond 
the  first  local  maxima.  Using  this  method  the  distance,  Rc.  becomes  6  cm.  In  both 
circumstances  the  distance  to  the  plume  is  beyond  the  value  of  Rc.  One  final  method  of 
evaluating  the  effectiveness  of  the  far  field  approximations  involves  looking  at  the 
directionality  of  the  transducer  at  the  range  of  interest  Using  Clay  and  Medwin23  the 
directionality  is  determined  and  plotted  for  the  distance  perpendicular  (rp)  to  the  propagation 
path: 

^  ^  2[j,(kasin<l))] 
kasin<() 

where  <|)  =arctan(rp/25cm). 


Perpendicular  range  from  center  of  plume  in  cm 
Figure  2.3 


22  Clarence  S.  Clay  and  Herman  Medwin,  Acoustical  Ocf  biography: 
Principles  and  Applications  (New  York:  John  Wiley  &  Sons, 1977),  p.l55. 
22  Clay,  p.  454. 
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This  indicates  that  for  the  plume  with  a  rough  diameter  less  than  1  cm  the  non  planer  terms 
are  small,  but  for  the  plumes  beyond  2  cm  in  diameter  the  value  falls  below  (1/2)'^.  This 
half  power  value  indicates  strong  non  planer  effects. 

The  signal  generation  and  data  acquisition  system  consisted  of  a  Real  Time  Systems 
(RTS)  model  4700/MR  Acoustic  Measurements  System,  using  an  arbitrary  wave  form 
generator  (AWG).  The  AWG  has  a  64K  memory  with  a  maximum  frequency  of  30  MHz, 
and  12  bit  amplitude  resolution.  The  repetition  rate  for  the  signals  generated  is  variable 
from  1  microsecond  up  to  10  seconds  in  1  microsecond  steps.  The  signal  acquisition 
system  uses  a  12  bit  digitizing  multiple  record  (MR)  board  with  a  maximum  sampling 
frequency  of  20  MHz  and  signal  length  of  up  to  1024  points  per  wave  form  and  store 
sequentially  up  to  2048  wave  forms. 

Temperature  measurements  of  the  temperature  profile  of  the  tank  during  the 
experiment  used  an  Omega  WB-AAI-B  High  Resolution  temperature  sampling  board  in  an 
NEC  Power  Mate  Portable  SX  with  8  channels  of  chromel-alumel  (K)  30  AWG 
thermocouple  wire  on  a  plastic  rod.  The  thermocouple  probes  were  spaced  approximately 
7.6  cm  apart  up  to  53.5  cm  with  the  one  probe  placed  at  the  1  cm  to  get  a  near  bottom 
sample.  Temperature  samples  were  taken  from  all  eight  channels  at  1  second  intervals. 
The  stated  accuracy  of  the  temporal  tc  board  was  0.7  °C  with  a  resolution  of  0.02  °C  for 
this  sample  rate;  but  all  thermocouples  were  verified  to  be  within  0. 1°  in  air  prior  to  the  start 
of  the  experiment. 

The  system  was  controlled  by  two  Macintosh  computers.  A  Ilfx  contained  the 
video  digitizing  board  and  was  used  mostly  to  obtain  the  video  data.  A  IIx  was 
programmed  to  run  the  RTS  system.  All  the  processing  of  the  data  from  this  facility  was 
done  on  a  Macintosh  IIvx. 
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2.2  Salt  Water  System 


The  other  laboratory  system  used  a  large  tank.  The  tank  had  cylindrical  ends  with 
radii  of  122  cm  with  a  122  cm  section  inserted  between  the  cylindrical  ends.  The  tank  had 
a  nominal  depth  of  1S2  cm.  This  gave  a  working  volume  of  1 1.65  m^.  An  anodized 
aluminum  frame  place  in  the  tank  held  a  heater  identical  to  the  one  used  in  the  other  lab.  In 
addition,  the  frame  held  rods  with  Teflon  transducer  holders  which  allowed  face  to  face 
transducer  distances  of  163  cm.  The  transducers  used  were  also  Panametrics  but  of  0.5" 
diameter.  This  diameter  gives  a  value  of  Rc=8.45  cm.  using  the  National  Standards 
Institute  method. 


Figure  2.4 

The  plot  of  the  directionality  demonstrates  that  even  for  plumes  greater  than  5  cm  lire  beam 
can  essentially  be  considered  a  plane  wave. 
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The  temperature  sampling  system  used  to  examine  the  nature  of  the  plume  consisted 
of  the  same  system  as  the  other  laboratory  except  that  the  array  consisted  of  4  groups  of  2 
thermocouples  spaced  10  cm  apart  and  placed  at  a  height  of  13  cm  above  the  heater  so  as 
not  to  interfere  with  the  signal  transmitted  at  5  cm.  The  system  did  not  give  an  exact 
qualitative  assessment  of  the  plume  at  the  height  of  acoustic  illumination.  It  did  give 
indications  of  the  trend  in  the  turbulence  and  provided  a  strong  indication  of  the  variation  in 
the  plume.  For  other  heights  above  the  plume  the  array  was  moved  up  as  appropriate. 

The  5  volt  peak  to  peak  995  kHz  signal  for  this  tank  was  generated  using  an 
Analogic  Polynomial  Waveform  Synthesizer  Model  2020.  The  signal  was  received  on  a 
Lecroy  9430  10  bit  150  MHz  digital  oscilloscope  with  a  memory  card  using  an 
external  clock  provided  by  an  HP  8 1 16A  signal  generator  set  at  10  MHz  to  allow  256 
repetitions  of  200  sample  points  to  be  stored  on  the  memory  disk  at  10  Hz.  The  data 
transmission  and  retrieval  system  was  controlled  by  a  H&H  Innovative  Solutions  486 
computer,  programmed  to  control  the  wave  form  synthesizer  and  digital  oscilloscope.  Data 
transmission  from  the  digital  oscilloscope  to  the  local  computer  network  was  accomplished 
using  a  program  called  ATF  developed  at  WHOI.  The  limit  on  the  storage  of  wave  forms  to 
256  was  overcome  by  running  the  experiment  twice  for  each  power  level  of  the  heater  to 
obtain  approximately  50  seconds  of  data. 
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Salt  Water 
Tank 


Figure  2.5 
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Chapter  3 


Theory 

The  experiment  described  in  the  previous  chapter  provides  an  easily  measured  result 
of  the  scattering  through  a  plume.  In  order  to  analyze  the  data,  several  aspects  of  the 
experiment  need  to  be  discussed.  The  nature  of  the  plume  and  the  experimental 
requirements  to  ensure  that  the  plume  can  be  properly  sampled  need  to  be  determined.  Once 
the  requirements  for  sampling  the  plume  are  established,  the  signal  to  be  used  must  be 
produced  to  achieve  the  desired  results. 

After  the  experiment,  the  analysis  of  the  data  requires  some  approximations.  The 
two  most  common  approximations  are  the  Bom  and  tl»  Rytov  approximations.  These  two 
approximations  are  explored  for  the  single  iteration  performance  of  each  to  assess 
applicability. 

Finally,  a  brief  explanation  of  the  scintillation  index  is  provided.  This  allows  a 
slightly  different  insight  into  the  effects  of  the  turbulence  on  a  signal,  and  may  assist  in 
visualizing  the  interaction  between  the  plume  and  the  signal. 

3.1  Buoyant  Plume 

In  order  to  properly  set  up  the  experiment,  the  rate  of  change  of  the  phenomena 
being  sampled  needs  to  be  determined.  The  actual  phenomena  need  not  be  understood  in 
order  to  view  it  Actually,  improper  sampling  of  the  phenomena  can  be  useful  if  it  is 
understood  what  the  sampling  limits  are.  The  plume  is  created  using  a  small  heat  source  on 
the  bottom  of  the  tank.  By  heating  the  water  around  the  source  a  narrow  buoyant  plume  is 
created.  The  nature  of  the  plume  is  determined  by  the  power  applied  to  the  heater.  At  the 
heaters  upper  power  limit  a  turbulent  plume  in  our  experimental  system  has  a  flow  rate  (U) 
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of  ®2.75  cm/s.  For  a  turbulent  system  the  model  for  the  turbulent  flow  provides  an  overall 
velocity  of  U+u,  where  u  is  a  random  variable  with  a  maxlul=aU  and  a<l.  The 
Kolmogorov  microscales  are  defined  for  length  (tj),  time  (t),  and  velocity  (\)).  As  such, 
the  minimum  inner  (small)  scale  for  a  turbulent  plume  in  water  is  given  by^'*  equation  (4); 
where  v  is  the  kinematic  viscosity  (1*10'^  m^s'l  for  water  at  20°C),  and  e  is  the 
dissipation  rate  per  unit  mass  (m^s'^).  The  expression  for  e  is  given  by^s  equation  (5). 


The  integral  (outer)  scale  of  the  plume  is  1,  and  u  is  the  random  variable  of  the  velocity.  By 
using  the  value  for  the  smallest  possible  eddy  size  and  the  velocity  of  the  plume,  with  the 
requirement  for  every  sample  to  be  no  more  than  one  eddy  size  from  the  previous  sample, 
an  effective  sample  rate  can  be  determined.  For  a  maximum  value  of  U+u=  5.5  cm/s  and 
1=1  cm:  h=88*10'^  m  and  a  velocity  of  5.5  cm/s  the  sample  rate  is  624  Hz  which 
corresponds  to  a  time  of  1.6  ms.  Initial  experimental  estimates  of  the  maximum  velocity  (2 
cm/s)  and  the  minimum  1  (2.5  cm)  gave  a  value  for  h=236*10‘6  m  and  a  value  for  the 
sample  rate  at  84.6  Hz.  The  experiment  was  designed  using  this  initial  value  of  a  sampling 
frequency.  Subsequent  experimentation  using  a  sample  rate  of  100  Hz  lead  to  the 
discovery  that  variations  in  the  plume  were  limited  to  less  than  2.5  Hz.  This  indicates  a 


H.  Tennekes  and  J.L.  Lumley,  A  First  Course  in  Turbulence  (Cambridge: 
The  MIT  Press  ,1972),  p.20. 

2^  Tennekes,  p.20. 
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value  of  u»2.3  mm/s  for  a  value  of  1  of  1.0  cm  and  corresponding  values  of  t|=  8  mm 
(using  the  same  method  as  above  and  2  cm),  x=.91  s  and  \)=1.1  mm/s . 

Nothing  else  is  necessary  to  establish  the  criteria  for  sampling  the  plume  created. 
Since  the  plume  is  initially  laminar  and  eventually  becomes  turbulent,  the  limiting  values  for 
the  experiment  are  outlined  above. 

The  structure  of  the  plume  itself  as  it  becomes  turbulent  is  of  interest  in  order  to 
investigate  the  trends  in  the  received  acoustical  data.  From  Batchelor^^  the  weight 
deficiency  produced  per  unit  time  at  the  source  (poF)  for  a  laminar  plume  is  given  by 
equation  (6). 


T(r)  -  T, 
To 


2Jtr  d  r 

j 


(6) 


Here  g  is  the  acceleration  due  to  gravity,  Tq  is  the  absolute  temperature  of  the  field,  and 
T(r)  is  the  absolute  temperature  of  the  plume.  T(r)  is  a  function  of  r;  the  axial  distance  from 
the  center  of  the  plume.  This  value  (F)  is  conserved  for  a  given  height  z.  A  Gaussian 
temperature  disuibution  is  given  by  equations  (7)  and  (8). 

il 

T(r)  -  To  =  ATe"’'  =  AT(r)  (7) 

F  =  (8, 

•^0 

Once  the  plume  becomes  turbulent  the  mean  temperature  distribution  continues  to 
remain  Gaussian.27  Batchelor  provides  equations  to  use  for  the  average  T(r)-To  and  mean 
velocity  ( U)  for  a  buoyant  plume  in  air  in  equations  (9)  and  (10). 


U  =  4.  7 


-1/3 


-96rVz 


(9) 


G.K.  Batchelor, “Heat  Convection  and  Bouyancy  Effects  in 
Fluids,  "£>uarteriy  Journal  of  the  Royal  Meterological  Society, SO  (1954), 
pp.  339-358. 

Batchelor,  p.344. 
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T(r)  -  T„ 


11  0  — 

gz^ 


(10) 


An  example  using  U=2  cm/s  and  a  AT  near  the  origin  of  10°  C  is  shown  in  figure  3. 1 : 


-0.1  -0.05  0  0.05  0.1 


Radius  from  center  in  meters 
Figure  3.1 

One  of  the  most  significant  features  of  the  plot  is  the  -5/3  exponential  decay  of  the 
peak  temperature  as  the  height  increases.  This  is  opposed  to  a  - 1  decay  noted  for  a  laminar 
plume.  This  decay  indicates  that  as  the  height  above  the  source  increases  the  scintillation 
index  will  decrease  at  a  lesser  rate  than  the  -5/3,  because,  the  peak  values  of  the  index  of 
refraction  will  decrease  at  a  lesser  rate  due  to  spreading  of  the  plume.  In  addition  the 
broadening  of  the  plume  indicates  a  broadening  of  the  eddies  and  demonstrates  that  the 
peak  modulation  frequency  of  the  received  signal  should  decrease  at  a  rate  corresponding  to 
a  broadening  rate  of  the  plume. 
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3.2  Signal  Description 

The  initial  genesis  of  the  problem  being  investigated  was  a  telemetry  question 
involving  the  possibility  of  discrete  turbulent  interference  of  shallow  water  signals.  The 
nature  of  this  interference  would  be  discrete  phenomena,  such  as  fresh  water  springs  or 
heat  vents.  The  common  effect  of  both  sources  of  interference  is  the  variation  in  the  index 
of  refraction  within  the  plume  by  variations  in  salinity  or  temperature.  Many  telemetry 
devices  involve  the  use  of  signals  in  the  lO’s  of  kHz  range.  An  initial  value  of  10  kHz  for 
a  signal  was  used  for  modeling.  The  natural  phenomena  discussed  above  are  in  general 
meter  scale  sources. 

A  common  method  for  evaluating  scattering  involves  the  factor  ka:  k  is  the  wave 
number  w/c  (w=2*7t*f,  c=  speed  of  sound  in  water)  and  'a'  is  an  effective  radius  of  the 
scatterer.  For  a  10  kHz  signal,  c=1500  m/s  and  an  ’a’=l  m  ka  becomes  =42;  which 
scales  the  frequency  to  1  MHz  for  this  heater  with  a  1  cm  plume. 

To  allow  easy  processing  at  the  desired  frequency  a  continuous  wave  (CW)  signal 
was  chosen.  The  length  of  the  signal  is  limited  by  the  necessary  sampling  rate  and  the 
avoidance  of  backscattered  echoes  off  of  the  tank  walls.  A  100  Hz  sampling  rate  limits  the 
total  length  of  the  signal  to  10ms.  To  avoid  echoes  in  the  reception  of  the  signal  the  signal 
length  should  be  less  than  the  transit  time;  the  transit  time  for  the  small  tank  is 
approximately  0.2  ms  (transducer  to  back  wall  and  back  to  transducer).  To  allow  a  large 
tolerance  factor  for  time  filtering  the  echo  signals  a  factor  of  20  reduces  the  limit  to  0.01  ms 
or  lOcycles  at  IMHz.  Ten  cycles  still  provides  a  strong  IMhz  signal  and  allows  about  10 
samples  per  cycle  with  a  sampling  rate  of  10  MHz. 

The  signal  in  Figure  3.2  is  the  raw  received  signal  for  the  experiment.  In  order  to 
determine  the  magnitude  and  phase  of  the  signal  most  accurately,  a  block  FFT  (i.e.  no 
windowing  except  the  basic  rectangular  window)  is  performed  and  the  result  is  shown  in 
Figure  3.3.  A  rectangular  window  provides  an  accurate  assessment  of  the  frequency 
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amplitude  with  one  frequency  present  with  a  processing  gain  of  l.fP* . 


Figure  3.2 


2®Fredrick  J.  Harris, “On  the  Use  of  Windows  for  Harmonic  Analysis  with 
the  Discrete  Fourier  Transform,*  Proceedings  of  the  IEEE,  66,  No . 1 
(January  1978), 54. 
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3.3  Interaction  of  Signal  and  Plume 

The  image  of  the  plume  given  by  the  acoustic  scattering  through  the  plume  is  really 
a  measure  of  a  slice  of  the  plume.  Since  the  signal  propagation  path  is  directly  through  the 
plume  the  resulting  signal  is  the  result  of  the  average  index  of  refraction  across  the  plume. 


3.3.1  Wave  Equation 

Starting  with  the  wave  equation  from  Ingard  and  Morse, and  following  the 
applicable  portions  of  the  derivation  in  Oeschger.^Oan  analytic  form  for  the  forward 
scattered  signal  is  obtained. 


c=  ae  ' 


(iyjr.  K(r,  t)V)  .  .  V  .  .  V  A)p(r,  t) 


p(  r ,  t )  is  the  pressure  at  location  r  and  time  t;  c  is  the  ambient  sound  speed; 

K  —  K 

=  - -  is  the  relative  compressibility;  is  the  ambient  compressibility; 

^0 

y  =  ^ ^  is  the  relative  density;  p,,  is  the  ambient  density;  and  u  is  the  fluid  velocity. 
Po 

Since  the  value  for  x  as  calculated  above  indicates  that  the  variability  of  the 
inhomogeneities  is  many  times  slower  than  the  period  of  the  sound  wave,  the  time 
dependence  will  be  ignored  for  Yk  and  Yp.  In  addition,  the  direction  of  u  is  essentially 
perpendicular  to  the  direction  of  the  sound  wave  propagation;  this  eliminates  the  last  two 
terms  in  the  parenthesis  of  (1 1).  Taking  the  Fourier  transform 


Philip  M.  Morse  and  K.  Uno  Ingard,  Theoretical  Acoustics  (Princeton, 
New  Jersey:  Princeton  University  Press, 1968) ,pp.  407-409. 

30  oeschger,pp.  12-15. 
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(p(r ,  w)  =  f  P(r ,  ty“"dt  and,  P(^/  t)  =  f  p(r ,  w)e’^“"dw) leaves 
equation  (12),  which  is  equation  8.1.12  from  Morse  and  Ingard.^^ 

(V  +  k^)p(r ,  w)  =  +  Y.  •  (Yp(r)Y))p(r ,  w)  (12) 

The  scattered  pressure  is  determined  by  the  integral  solution  to  (12),  which,  for  the 
boundary  free  case,  is  given  in  Morse  and  Ingard  as  equation  8.1.13^2 

pjr  ,  w)  =  (13) 

J  {k^Y^(r  ')p(r  '  ,  w)g(r|  r  ■)  +  Yp(j^  ')Yp(r  '  ,  w)  •  Vg( r  I  r  ■)}dr  ' 

The  relative  compressibility  and  the  relative  density  are  the  two  primary  sources  of 
variation  in  the  effective  index  of  refraction  in  fluid  plume.  The  relative  importance  of  these 
two  factors  can  be  determined  by  considering  the  variation  in  temperature  in  the  plume. 
Using  a  value  of  5  °C  for  the  change  across  the  plume;  the  values  for  the  change  in  density 
and  sound  speed  for  this  variation  can  be  found.^^  The  sound  speed  at  20  °C  is  1482.3 
and  at  25  °C  it  is  1497.5.  The  value  for  the  density  of  fresh  water  at  20  °C  is  998.2 13  and 
at  25  °C  is  997.053. 

2  2 

Y^  =  ~  ^‘^  ■=■0.019  (14) 

pc 

y.  =  P  =-0.00 12 

P 

The  ratio  between  y  k  and  Yp  is  about  15.  The  scattering  due  to  changes  in  the 
compressibility  is  far  more  significant  than  that  due  to  variations  in  the  relative  density. 

^^Morse,  p.410. 

^^Morse,  p.410. 

L.E.  Kinsler  et  al Fundamentals  of  Acoustics  (New  York:  Wiley  and 
Sons, 1982 ) ,  p .  107 . 
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An  expression  for  in  terms  of  temperature  in  accordance  with  Oeschger^  ^  can  be 


found  by  assuming  that  p^^po  ^d  using  a  linear  approximation  of  the  sound  speed  such 
dc 

that  c  =  — At  +  Cq.  Since  c«co,  the  expression  for  Yk  becomes  equation  (15)which 
dT 

reduces  to  equation  (16). 


Yk  = 


(15) 


Yk 


2  dc  .  . 
- —  AT(r). 


Co  dT 


(16) 


Using  these  approximations,  the  scattered  equation  becomes  equation  (17),  where 
dc 

dT’ 


Pg(r,w)  =  -2ak^  J  AT(r ')p(r  '  ,  w)g(r I  r  Odt' '  (17) 


The  generalized  asymptotic  form  of  the  Greens  function  for  one,  two  and  three 
dimensions  from  Oeschger^s  reduces  to  equation  (18)  where  d  is  the  dimensionality,  and 
I  r  I  is  the  radial  coordinate  in  three  and  two  dimensions  or  the  linear  coordinate  of  the  one¬ 
dimensional  case. 


^  .fill)  f±±)  I 
(r  I  r  ')  =  (nl  r  |)  2  2  ^  ^  ^  ^  exp 


i  kl  r  I  -k. 


•  r  ' 


(18) 


^^oeschger,  p.  14. 
Oeschger,  p.l5. 


31 


Figure  3.4 


3.3.2  Born  Approximation 

In  order  to  obtain  a  solution,  a  single  iteration  Bom  approximation  will  be  used.  In 
accordance  with  Ishimaru^^,  in  order  to  use  the  Bom  approximation  successfully:  (n^- 
l)kD«l.  This  is  satisfied  for  a  plume  up  to  1.2  cm  in  diameter.  The  Bom  approximation 
allows  the  substitution  inside  the  integral  for  p*  of  pi(w)exp(i  •  r  ’)  for  p(  r  ',w). 

Upon  making  this  substitution  the  ,'  cattered  field  becomes  equation  (19). 

P3(r,w)  =  a(2Jir)Tk^^^p.  (w)  exp|^i|^kr  + 

Xj  AT(r  ' 

The  integral  portion,  for  ks=ki,  reduces  to  J  AT(r  ')dr  Equation  (20)  is  in  cylindrical 
coordinates  for  a  two  dimensional  system  withAT(r')  =  ATe 

Ishimaru, p . 22 . 
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Pg{r  ,  vi)  =  oATb^  (w)  exp^i^kr  +  j  ^20) 


Equation  (20)  in  terms  of  the  total  field  Pt=pd+Ps  is  shown  in  figure  3.5  where  pd  is  the 
direct  unscattered  field. 


3.3.3  Rytov  Approximation 

Some^'^  recommend  the  Rytov  method  for  solving  in  the  forward  scattered 
direction.This  method  involves  the  substitution  of  p(r,w)=pi(w)exp(i\|i(r));  where 
\|/(r)=S(r)-i*log(p(r)/pi).  Equation  (21)  is  equation  (12)  for  the  case  of  Yk  >>  Yp- 

(V^  +  k^)p(r,w)  =  -k^y^{x)p{x  ,vi)  (21) 


Ishimaru,  p.  349. 
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Substituting  p(  r  ',w)=pi(w)exp(i\/(  r  '))  and  then  following  the  substitutions  for  the 
derivation  of  the  Rytov  approximation  in  Ishimani^®  a  solution  results.  Using 
V  =  Vo  +  ¥i  with  vitg  from  =0  the  form  becomes  equation  (22). 

(V  +  k')p,(r,  w)v,  =  -(VxjTj  •  V\|/,  +  k^yjpjr,  w)  .  (22) 

Equation  (22)  is  in  the  form  of  the  inhomogeneous  wave  equation  which  when  converted  to 
integral  form  for  becomes  equation  (23). 

(23) 

Vi(r,  w)  =  — - - -  f  (VVj  ■  VVi  +  k'Yjpo(r  '  ,  w)g(r|  r  ')dr  ' 

Equation  (23)  yields  equation  (24)  for  a  first  iteration  of  Vj  =0  on  the  right  hand  side. 
Vi(r,w)  =  — - -  f  k^YKPo(*^  ’  /  w)g(r|  r ')dr  '  (24) 

Since  po  is  defined  for  Yk=^  ^T(  r )  =0  for  all  r,  the  variations  in  Po  are  due  only  to 
propagation  in  the  media.  The  distance  across  the  plume  is  only  about  2.5  cm  and  the 
decay  from  r=25  cm  to  r=50  cm  in  the  signal  is  only  4%.  Using  p^  for  the  magnitude  of 
Po  ( r ,  w )  and  Po  ( r  '  ,  w )  is  not  a  bad  approximation. 

VoCr)  =  S(r)  -  ilog(p(r)/p.  ) 

S(r)=k,  •  r  -  ko 

Vio(i^)  =  e'^'  ''g(r|  r ')k^YKdr  '  (25) 

For  Yk  =  -2  aAT(  r  ' )  and  the  Green's  function  as  defined  in  (18),  equation  (25) 
arrives  at  equation  (26). 


^^Ishimaru,  p.  350. 


d  +  l 


Vio  (r  ')  =  oi(2Jt|r'|)2k^e^*'^e'  ‘  * 


Equation  (27)  is  equation  (26)  in  cylindrical  coordinates. 


I  AT(r  ••)dr  “  (26) 


\|/io(r')  =  txATb' 


21  r'  I 


))• 


(27) 


Using  AT(r  ')  =  ATe  the  result  becomes  equation  (28). 


V„(r')  - 


(28) 


Directly  inserting  this  into  ps(r)=pi(w)e't'‘“  gives  a  similar  solution  to  the  Bom 
approximation  but  with  slightly  lower  slope  and  an  opposite  phase  variation.  The  total 
received  signal  pt  is  shown  in  figure  3.6. 

If  the  asymptotic  solution  for  the  Green’s  function  given  in  (15)  cannot  be  used  for 
the  determination  of  \|/io(r),  then  for  the  axisymmetric  case  the  Green’s  function  is  given 
by  equation  (29),  and  the  expression  becomes  equation  (30). 

g(r|  r  ■)  =  i  H‘'’(k|r  -  r 'I).  (29) 

4 

\|/j<,(r)  =  'I  AT(r ')e‘*‘  '' H<'’(k|r  -  r 'Ddr ’.  (30) 

This  function  has  a  strong  pole  at  the  points  where  r=r’.  To  take  the  Rytov  approximation 
to  the  next  approximation  beyond  Vio  to  Viio  insert  Vio  as  shown  in  equation  (31). 

(31) 

Viio(r)  =  +  —r^ — :  1  G(r  -  r  ’)V\jrjo(r  ')  •  V\|ri(,(r  ')p(r  '  ,  v/)dr  ’ 

p(r ,  w)  •' 
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Figure  3.6 


3.4  Cylindrical  Approximation 

In  order  to  make  the  comparison  between  different  power  levels  of  the  plume,  the 
plume  will  be  modeled  as  a  fluid  cylinder  of  constant  density.  From  Sianion^^  the  value 
for  the  scattered  pressure  field  is  given  by  equations  (32)  and  (33). 

=  Po  i""  cos  (m(|))j'„(kr) 

m  sO 

Pine  =  Po  i  COS (m<)))j„(k’r)  (32) 

m*o 

Ps  =  Po  i  cos(m(t>)  H<‘^kr) 

m®0 

T.  K.  Stanton, "Sound  Scattering  by  Cylinders  of  Finite  Length. I. 
Fluid  Cylinders, " The  Journal  of  the  Acoustical  Society  of  America,  83 , 
No.l  (January  1988),  55-63. 
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Phase  Magnitude 


The  Neumann  numbers  are  Eo=i  =1  and  Cm  =2  for  all  other  m;  and  k*  is  the  wave  number 
inside  the  cylinder. 

Bm=-eni/'"/(l+/Cn,) 


W„(ka)]y 

- 

'  (ka)/J„  '  (ka)] 

-'’m  '  (k*a 

)  J™(ka)]/[j„(k*a; 

'(ka) 

-  gh 

p. 


The  response  of  the  scattered  field  in  the  forward  direction  for  the  cylinder  model  divided 
by  Po  is  shown  in  figure  3.7. 


Figure  3.7 
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The  phase  response  of  both  the  cylinder  model  and  the  Bom  approximation 
decreases  with  increasing  AT,  but  the  phase  of  the  Rylov  approximation  increases  with 
increasing  AT . 

The  cylinder  can  be  made  to  match  any  of  the  given  magnitudes  for  different 
temperatures.  For  the  Bom  approximation  the  ratio  is  approximately  8. 16*Cyl/Bom.  The 
cylinder  is  used  only  to  allow  a  common  base  for  comparison  of  the  temperature  profiles  in 
each  experiment  This  common  base  allows  for  analysis  of  the  mean  temperature  variations 
within  the  plume. 

3.5  The  Scintillation  Index 

The  Scintillation  index  is  a  measure  of  the  variation  of  the  plume. 


In  this  circumstance  the  intensity  is  going  to  vary  with  the  variation  in  the  index  of 
refraction.  For  low  powers  the  plume  varies  very  little  and  as  such  the  0^=0  (or  very 
small).  As  the  variation  in  the  index  of  refraction  varies  in  magnitude  more  and  more  as  the 
turbulence  increases  the  increases  as  well.  At  some  point,  if  the  plume  becomes  fully 
turbulent,  the  will  approach  a  constant  that  will  be  a  measure  of  the  maximum  variation 
in  the  index  of  refraction  achieved  in  plume.  It  may  also  show  a  decrease  in  as  power  to 
the  heater  continues  to  increases  due  to  the  lowering  of  the  point  at  which  the  plume 
becomes  turbulent  and  the  additional  height  with  a  (-5/3)  exponential  temperature  decrease. 
This  may  be  overcome  by  the  effects  of  spreading.  Looking  at  the  plume  as  a  divergent 
lens  the  normalized  magnitude  should  not  return  to  unity  with  increasing  turbulence. 

In  this  model  of  a  turbulent  plume  if  the  peak  increase  in  temperature  in  the  plume  is 
near  10°C  and  there  is  assumed  to  be  a  Gaussian  distribution  in  lime  with  no  spreading,  the 
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will  approach  SXIO'^.  This  value  is  based  on  the  Rytov  magnitude  and  that  the 
measured  peak  temperature  is  approximately  5X  the  mean  difference  from  the  zero  power 
value  of  temperature. 

The  scintillation  index  provides  a  little  more  insight  into  the  turbulent  plume.  If 
there  is  no  variation  in  the  plume  the  scintillation  index  will  be  zero.  However,  as  the 
variance  increases  the  scintillation  index  increases.  In  our  experiment,  with  weak 
scattering,  the  amplitude  modulation  of  the  signal  is  small,  so  the  scintillation  index  can  be 
seen  as  a  measure  of  the  mean  modulation  index  of  the  signal.  This  analogy  allows  the 
visualization  of  the  signal  as  a  carrier  signal  with  a  modulated  signal  produced  by  the 
variations  in  the  index  of  refraction  and  added  on  to  the  carrier. 
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Chapter  4 


Experimental  Results 

The  results  of  the  previous  section  were  used  to  look  at  the  forward  scattering 
through  the  thermal  plume  generated  in  the  experiments  also  previously  described.  The 
results  from  the  freshwater  experiment  are  compared  to  the  expected  results,  and  then  they 
are  extended  to  a  more  complicated  experiment  Each  of  the  different  heights  were 
conducted  at  different  times  and  under  slightly  different  tank  conditions.  The  15  cm 
experiment  was  conducted  at  the  same  height  as  the  plume  temperature  measurements. 
Because  of  this,  the  15  cm  data  will  be  used  for  specific  comparisons  to  the  theoretical 
models.  Extensions  to  the  other  heights  will  concern  the  nature  of  the  plume. 

4.1  Fresh  Water  Tank 
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Temperature  in 


15.2  cm  Temperature  response 


Figure  4.2 

The  temperature  profile  for  the  tank  is  shown  in  figure  4. 1  for  the  duration  of  the 
experiment  Using  a  10  second  average  the  temperature  difference  over  the  duration  of  the 
experiment  is  0. 107 1°  C.  A  plot  of  the  temperature  change  of  the  thermocouple  at  1 5.2  cm 
is  shown  in  figure  4.2. 


4.1.1  Magnitude  and  Phase  Mean  Response 

The  temperature  data  at  60  volts  (1 1.5  Watts)  in  Oeschger  v  a  taken  at «  15  cm. 
The  variation  in  temperature  with  power  was  found  to  be  linear  with  power^®.  Figure  4.3 
is  a  plot  of  the  magnitude  response  of  the  Rytov  method.  The  Bom  approximation,  and  the 
results  for  the  signal  at  15  cm.  The  negative  signs  indicate  decreasing  power  series. 


^^Personnel  conversation  with  Oeschger,  24  June  1993. 
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Figure  4.3 

The  amplitude  response  seems  to  indicate  a  lag  time  in  the  plume.  The  plumes  were 
allowed  to  stabilize  5  minutes  prior  to  each  signal  set  and  were  visually  checked  for 
consistency  when  the  plume  became  visible.  The  magnitude  response  of  the  plume 
parallels  the  Rytov  approximation  initially  but  more  closely  approaches  the  Bom  as  the 
power  is  decreased.  The  response  follows  a  similar  pattern  as  the  power  is  increased.  A 
non-linear  temperature  vs.  power  response  at  low  powers  could  indicate  a  higher 
temperature  than  the  linear  power  vs.  temperature  model  provides.  The  higher  powers  are 
also  subject  to  increasing  turbulence  which  changes  some  of  the  heater  energy  into 
turbulent  motion  and  accounts  for  some  of  the  flattening  of  the  amplitude. 

An  effort  to  correct  for  the  phase  changes  due  to  temperature  variations,  needs  to 
remove  the  phase  due  to  increased  travel  time  from  the  0.1071°C  increase  in  ambient 
temperature.  The  slope  of  the  phase  change  with  power  of  the  heater  is  equal  to  .0014. 

This  was  obtained  using  the  mean  slope  for  the  decreasing  power  series  and  the  mean  slope 
for  the  increasing  power  series  and  removing  the  average  increase  to  obtain  the  slope  due  to 
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the  changing  power  level  alone. 


Figure  4.4 

Although  it  appears  that  the  Bom  approximation  provides  a  better  match  of  the 
amplitude  response  of  the  signal  than  the  Rytov  approximation,  the  phase  response  is 
clearly  better  matched  by  the  Rytov  approximation. 

Figure  4.5  shows  the  magnitude  responses  for  each  of  the  heights  used.  If  the 
magnitude  responses  for  the  different  heights  are  compared,  the  effects  of  the  decay  in 
temperature  with  the  height  of  the  plume  can  be  seen  as  in  figure  4.5.  Figure  4.6  shows 
this  variation  in  terms  of  temperature  of  the  cylinder. 

According  to  discussion  in  chapter  3  the  peak  temperature  should  decrease  linearly 
with  height  for  a  laminar  plume.  Even  though  the  heights  used  were  done  at  different  times, 
the  laminar  plume  equivalent  temperature  model  plotted  as  a  function  of  height  for  the  lower 
po  .vers  (before  the  30  cm  signal  becomes  turbulent),  as  in  figure  4.7,  shows  temperature 
decay  as  nearly  linear  with  a  mean  slope  of  -.0236°Cycm. . 
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Cylinder 


Figure  4.6 
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For  the  turbulent  case,  which  is  most  nearly  achieved  only  for  the  highest  power 
level  based  on  video  and  the  scintillation  index,  a  similar  comparison  can  be  determined. 
The  radius  is  supposed  to  increase  linearly  with  height,  as  such,  the  variance  should  spread 
linearly.  The  temperature  decay  for  a  turbulent  system  is  nearly  -5/3  in  height,  a  net  -2/3 
magnitude  vs.  height  exponential  decay  in  signal  response  is  expected.  However,  the 
limited  number  of  available  turbulent  signal  responses  ensures  a  skewed  result.  To  match 
the  signal  at  5  cm  a  factor  of  0.2535  is  needed.  Figure  4.8  shows  the  plot  of  the  highest 
power  responses  for  each  of  the  heights.  The  signal  at  15  cm  appears  not  to  have  been  fully 
turbulent  at  the  power  level  shown.  The  error  at  30  cm  will  be  significant  for  the  constant 
used  at  5  cm.  Clearly  there  is  not  enough  data  to  do  a  good  comparison. 
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4.1.2  Frequency  Response 

The  propagated  signal  time  series  with  the  signal  sampled  at  100  Hz  when  the 
plume  is  transitioning  to  turbulent  appears  in  figure  4.9.  A  Fourier  transform  of  the  signal 
in  figure  4.9  illustrates  that  most  of  the  signal  is  well  below  2  Hz  and  is  shown  in  figure 
4.10.  The  sampling  rate  was  subsequently  reduced  to  10  Hz  which  allowed  a  closer  look 
at  the  low  frequency  modulation  of  the  signal.  This  sample  rate  (actually  the  pulse 
repetition  rate  of  the  signal  used)  lowered  the  highest  discernible  frequency  to  5  Hz  and 
reduced  the  frequency  bin  size  from  0.098  Hz  per  bin  to  0.0098  Hz  per  bin. 

If  the  FFT  of  the  signal  magnitude  squared  at  30  cm  is  examined,  and  a  rough  x2  fit 
is  made,  a  close  approximation  to  the  x2  distribution  can  be  seen.  Figure  4. 1 1  uses  a  ^2 
distribution  of  order  two  which  indicates  an  underlying  Gaussian  distribution  of  the  real 
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Normalized  frequency  response  Normalized  Magnitude  Squared 


and  imaginary  parts  of  the  received  amplitude  signal. 


Figure  4.9 
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The  x2  distribution  arises  from  the  sum  of  the  squares  of  the  two  Gaussian  variables  and 
is  of  the  form  of  an  exponential  decay.  The  lower  frequencies  dominate  this  spectrum  of 
signal  response.  Using  0.5  Hz  as  a  rough  estimate  of  the  size  of  the  inner  scale  at  0.707 
Hz  is  3.8  cm.  This  indicates  that  the  larger  eddies  dominate  the  modulation  of  the  signal. 

If  the  of  the  signal  for  several  heater  powers  is  shown,  the  increasing  modulation 
as  the  turbulence  increases  can  be  seen  as  in  figure  4.12.  The  signal  at  15  cm  never  fully 
became  turbulent  as  shown  in  the  plot  of  the  temperature  vs.  height  plot  in  figure  4.8  and  as 
will  be  seen  in  the  discussion  of  the  scintillation  index.  The  figure  does  demonstrate  the 
initial  appearance  of  higher  frequencies  before  the  low  frequencies  become  dominant. 
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Power  Level  of  Heater  in  Watts 


Propagated  Signal  in  Frequency  shifted  by  heater  power 


4.1.3  Scintillation  Index 

The  scintillation  index  was  determined  for  each  power  level  for  each  of  the  heights. 
Figure  4. 13  illustrates  the  signiflcant  points  discussed  in  chapter  3.  As  the  plume  becomes 
turbulent,  departing  from  the  laminar  plume,  the  scintillation  index  starts  to  increase.  The 
scintillation  inctex  increases  until  it  reaches  a  plateau,  where  the  turbulence  is  fully 
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developed.  In  addition  the  plume  may  not  be  turbulent  at  all  heights  above  the  source  but 
may  become  turbulent  at  some  point  above  it 

The  5  cm  scintillation  index  indicates  that  the  fluid  in  the  tank  was  not  fully  stable. 
The  scintillation  index  reaches  a  small  plateau  for  power  equal  to  0  watts  indicating  some 
turbulence  present  either  in  the  motion  of  the  water  in  the  tank  or  due  to  some  disruption  of 
a  stratification  of  the  tank.  However,  this  disruption  was  not  as  significant  as  the  plume 
turbulence  generated  at  higher  powers.  This  slight  turbulence  does  lead  to  the  similar  point 
of  departure  from  laminar  in  the  5  cm  plume  and  in  the  IS  cm  plume.  This  disruption  also 
leads  to  the  leveling  of  the  scintillation  index  at  a  lower  level  than  would  be  expected  due  to 
the  higher  peak  temperatures  at  5  cm  than  at  30  cm. 


xlO-3  Scintillation  Index  vs  Power 


Figure  4.13 

The  15  cm  scintillation  index  never  reaches  a  plateau  and  thus  the  15  cm  plume 
never  reaches  full  turbulence.  The  scintillation  index  does  indicate  that  some  sort  of 
disturbance  has  been  developed  as  the  experiment  has  proceeded  as  shown  by  the  lower 
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slope  for  the  increasing  power  portion  of  the  experiment.  The  IS  cm  scintillation  index  also 
indicates  that  the  higher  peak  temperature  within  the  plume  than  at  30  cm  leads  to  a  higher 
scintillation  index.  The  peak  scintillation  index  of  SxlO-^  corresponds  to  a  peak  of  25  °C. 
This  is  a  fairly  high  value  and  is  based  on  the  Rytov  approximation  and  the  5x  the 
equivalent  temperature  for  the  mean  difference  from  the  no  power  case  to  obtain  the  peak 
temperature.  The  Bom  approximation  would  lead  to  a  peak  temperature  of  10  °C.  That 
creates  a  reasonable  range  for  estimation  of  the  peak  temperature  from  the  scintillation 
index.  If  the  variance  of  the  plume  thennal  profile  increases  by  a  factor  of  9  (b  increases  by 
a  factor  of  3,  physically  spreads  by  a  factor  of  3)  the  temperatures  become  0.4  °C  for  the 
Bom  and  2.8  °C  for  the  Rytov. 

Since  the  peak  temperature  at  IS  cm  has  been  measured  and  is  at  least  S.2  °C  at  1  l.S 
Watts,  the  factor  by  which  b  has  increased  for  this  case  is  at  most  2.2.  for  the  Rytov 
approximation  and  1.3  for  the  Bom  approximation.  Unfortunately,  precise  statements 
caxmot  be  made  since  the  experiment  did  not  include  a  detailed  measurement  of  the  plume 
temperature  distribution  at  all  power  levels.  The  variations  in  the  scintillation  index  are  due 
in  general  to  either  the  increased  radius  of  the  plume  or  to  an  increase  in  peak  temperature 
of  the  plume. 

The  scintillation  index  for  30  cm  not  only  reaches  the  plateau  it  even  displays  a 
possible  decrease  in  the  scintillation  index  due  to  a  lower  height  for  the  point  of  departure 
from  laminar.  This  decrease  is  briefly  discussed  in  the  previous  chapter. 

The  freshwater  tank  data  provides  an  a  illustration  of  the  plume  in  both  its  laminar 
and  turbulent  states.  The  general  trends  in  the  effect  the  changing  state  of  the  plume  has  on 
the  propagated  signal  are  discussed  and  illustrated.  This  information  is  of  little  practical 
value  if  they  cannot  be  used  in  a  more  realistic  application. 
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4.2  Comparison  to  Salt  Water  Tank  Data 


4.2.1  Magnitude  and  Phase 

The  above  discussion  of  the  plume  in  the  fresh  water  tank  provides  some  guidance 
for  the  anticipated  results  for  a  thermal  plume.  The  small  fresh  water  tank  provides  a  well 
controlled  medium  in  which  tc  create  a  controlled  plume.  In  the  real  world  outside  of  this 
tank,  things  are  not  and  cannot  be  controlled.  As  such,  the  factors  by  which  we  have 
estimated  the  parameters  of  the  plume  may  not  provide  the  same  insight  as  in  the  controlled 
situation. 

The  large  tank  used  in  this  experiment  is  subject  to  stratification  in  salinity  and  in 
temperature.  In  addition  the  orientation  of  the  tank  and  the  presence  of  rir  conditioning  at 
one  end  of  the  tank  creates  internal  waves  within  the  tank.  These  internal  waves  greatly 
affect  the  propagated  signal  over  time.  The  tank  also  has  some  suspended  particulates  and 
gas  bubbles  which  can  greatly  affect  the  signal. 

The  effects  of  the  above  phenomena  can  be  seen  in  the  signals  and  are  easily  seen  in 
the  transmitted  pulses.  The  particulates  create  discrete  drop  outs  in  the  time  series 
magnitude  of  the  transmitted  signal.  Both  the  particulate  and  the  bubbles  have  such  sharp 
discontinuities  in  the  coefficient  of  refraction  that  the  propagated  signal  is  greatly 
diminished.  The  internal  waves  within  the  tank  are  of  a  longer  time  scale  than  the  sample 
time  periods  used  in  the  experiment.  An  example  of  the  received  magnitude  for  a  single 
pulse  train  is  shown  in  figure  4.14  for  a  44  watt  l^ter  power  and  a  height  of  15  cm  above 
the  plume. 

In  the  experiment  conducted  in  the  large  tank  the  healer  power  was  initially  at  zero 
watts  and  then  increased  to  the  peak  value  near  60  watts.  Then  the  power  was  reduced. 
This  is  the  opposite  of  the  method  used  in  the  smaller  tank.  The  plots  are  reversed  to 
present  the  same  orientation  as  the  previous  plots.  So  that  for  the  decreasing  power  series 
the  power  has  a  negative  sign  in  front. 
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Normalied 


Figure  4.14 


The  magnitude  response  of  the  signal  at  5  cm  in  height  are  shown  below  in  figure 


4.15. 


Figure  4.15 

This  plot  is  shown  first  to  illustrate  some  of  the  problems  inherent  in  the  use  of  a 
more  realistic  media.  The  large  magnitude  decrease  is  probably  due  to  the  passing  of  a 
large  eddy  between  the  source  and  the  receiver.  Unfortunately,  it  didn't  return.  The  plots 
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organized  in  the  same  manner  as  previous  plots  are  shown  below. 


Figure  4.16 


Power  in  Watts 


Figure  4.17 

The  phase  plot  does  not  show  the  effect  of  the  eddy,  but  it  does  not  show  the  effect  of  the 
heater  raising  water  temperature  either.  The  phase  seems  nearly  random  but  closer  to 
constant  at  the  lower  powers  as  it  should  be.  The  higher  power  phases  appear  to  be  much 
more  variable.  The  scale  of  the  variability  in  phase  is  more  than  ten  times  that  of  the  small 
tank. 

The  magnitude  plot  for  the  IS  cm  experiment  is  plotted  along  with  the  Bom  and 
Rytov  approximations  in  figure  4.18.  The  data  does  not  appear  to  have  been  subject  to  an 
eddy  but  merely  to  the  smaller  scale  fluctuations  in  the  media. 
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1 


Magnitude  plot  of  IScm  data  with  Bom  and  Rytov 


Power  in  Watts 
Figure  4.18 


It  appears  that  the  Bom  and  Rytov  approximations  are  still  valid  when  the  system  can  have 
some  of  the  larger  disturbances  removed.  The  phase  however  continues  to  be  nearly 
random  as  shown  in  figure  4. 19.  The  phase  varies  over  the  entire  27t  range  with  very 
small  variations  in  power. 
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The  plot  of  magnitude  at  30  cm  in  figure  4.20,  allows  some  of  the  decrease  in  amplitude 
with  increasing  power  to  be  seen,  but  also  shows  the  effect  of  becoming  turbulent  at  a 


lower  power  than  the  15  cm  case. 


Phase  Plot  of  30  cm  data 


The  phase  continues  to  be  diflicult  to  resolve.  The  phase  response  at  30  cm  is 
shown  in  figure  4.21.  The  unwrapped  phase  is  shown  in  Hgure  4.22  for  the  30  cm  case 
and  demonstrates  the  range  over  which  the  phase  changed  throughout  the  experiment  This 
plot  seems  to  illustrate  the  similar  increase  in  phase  as  the  experiment  progressed  as  seen  in 
the  smaller  tank.  However,  the  increase  in  phase  for  the  small  tank  was  only  0.8  radians 
during  the  15  cm  experiment  This  increase  in  phase  is  over  4  radians  and  due  to  daily 
heating  in  the  room.  In  order  to  remove  this  trend,  the  temperature  of  the  entire  tank  would 
have  to  be  continuously  measured  and  carefully  removed  from  the  trend  and  variations  in 
sound  speed  due  to  variations  in  the  heating  rate.  The  shear  magnitude  of  the  phase  change 
and  the  accuracy  of  the  measurements  over  this  range  would  have  to  be  fine  enough  to  see  a 
0.1  radian  phase  variation. 

4.2.2  Frequency  Response 

Another  measure  of  turbulence  is  the  modulation  of  the  magnitude  time  series.  In 
the  small  tank  a  x2  distribution  was  used.  The  x2  distribution  is  shown  with  the 
normalized  FFT  of  the  magnitude  squared  of  the  signal  from  the  salt  water  tank  in  figure 
4.23.  Figure  4.23  shows  a  x2  distribution  and  the  dominance  of  the  lower  frequencies 
within  the  spectrum.  The  fit  is  not  nearly  as  good  as  for  the  fresh  water  tank,  but  this  can 
mosdy  be  attributed  to  the  limited  turbulent  data  at  15  cm  in  the  salt  water  tank. 

A  Fourier  traasform  of  the  phase  response  helps  illustrate  the  lack  of  obvious  phase 
modulation.  As  shown  in  the  magnitude  and  phase  plots  the  phase  does  not  show  as 
simple  a  response  as  the  magnitude.  The  frequency  plots  of  the  magnitude,  figure  4.23, 
and  the  phase,  figure  4.24,  clearly  demonstrate  this  fact. 
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The  development  of  the  turbulence  can  be  seen  fairly  clearly  in  the  overlaid  plots  of 
frequency  distribution  vs.  power  of  the  magnitude  response  in  Hgure  4.25.  The  higher 
frequencies  tend  to  form  first  and  then  the  lower  frequencies  start  to  dominate  the  spectrum. 
This  indicates  that  the  plume  appears  to  break  up  initially  with  smaller  eddies  breaking  off 
first.  Once  the  plume  became  mote  turbulent  the  larger  eddies  form  and  dominate  the 
spectrum.  When  the  plume  is  fully  turbulent  the  spectrum  looks  to  be  closely  Gaussian. 

Fourier  Transform  of  Magnitude  Response 


Figure  4.25 
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The  frequency  spectrum  of  the  modulated  magnitude  time  series  appears  to  be  an 


excellent  indication  of  the  state  of  the  plume.  The  development  of  the  plume  as  power 
increases  can  also  be  seen  from  the  spectrum.  This  method  of  analysis  appears  to  be  a 
fairly  robust  measure  of  the  nature  of  the  plume. 


4.2,3  Scintillation  Index 

The  scintillation  index  was  also  determined  in  the  salt  water  tank  for  each  height. 
The  scintillation  index  for  similar  heights  for  the  two  tanks  are  overlaid  on  the  plots  below. 


Run  number  (power  1=0  watts,  power  15=57  watts,  power  29=0  watts) 

Figure  4.26 


The  scintillation  index  for  the  small  tank  is  generally  larger  as  expected  due  to  the 
increased  amount  of  disturbances  in  the  larger  tank.  However,  the  values  are  very  close 
between  the  two  tanks.  Figure  4.26  shows  the  5  cm  scintillation  response. 
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Run  number  (power  1=0  watts,  power  15=57  watts,  power  29=0  watts) 


Figure  4.27 


1 0-3  Scin  tillation  Index  at  30cm 


Run  number  (power  1=0  watts,  power  15=57  watts,  power  29=0  watts) 


Figure  4.28 


Figure  4.27  confirms  our  earlier  discussion  of  the  IS  cm  plume  in  the  small  tank.  It 
appears  that  the  plume  in  the  small  tank  never  reached  a  fully  turbulent  state.  From  the 
figure  4.27  the  plume  in  the  large  tank  may  also  not  be  fully  turbulent 

The  30  cm  plot  in  figure  4.28  illustrates  the  essentially  saturated  turbulence  for  the 
higher  power  levels.  The  plot  also  demonstrates  the  slightly  reduced  values  for  the  larger 
tank. 

Overall  the  scintillation  index  represents  one  of  the  most  robust  and  easily  achieved 
measurements  of  the  plume  structure  discussed  in  this  p^r.  Although  the  magnitude  of 
the  received  signal  can  be  reconstructed  with  careful  analysis  of  the  environmental  data,  it 
requires  extensive  manipulation.  It  also  appears  very  difficult  if  not  impossible  without 
highly  accurate  measurements  of  the  environment  to  recover  the  phase  relationships.  The 
frequency  spectrum  of  the  modulation  appears  to  be  a  good  indication  of  the  nature  of  the 
plume,  and  is  similar  to  the  scintillation  index  in  that  it  makes  use  of  the  variations  in  the 
time  series  of  the  magnitude. 

4.3  Scattering  Function 

In  order  to  use  the  information  within  this  thesis  for  telemetry  problems  the 
scattering  function  of  the  plume  needs  to  be  determined.  The  scattering  function  is  usually 
defined  in  terms  of  s(o),r)  and  can  be  considered  the  average  scatterer  cross  section^  ^  for 
each  pair  of  frequency  and  ranges.  In  this  case  only  the  frequency  distribution  was 
measured  and  as  such  the  scattering  function  will  be  shown  as  s(©).  From  Kennedy^ 2 
s((i))  can  also  be  seen  as  the  frequency  distribution  of  the  correlation  function,  R(t,x),  of  the 
received  signal.  This  is  sometimes  termed  the  power  spectral  density  S(o)).  The  mean 
correlation  function  of  the  signals  from  the  fresh  water  tank  at  30  cm  are  shown  in  figure 
4.29.  The  power  spectrum  is  also  defined  as  the  Fourier  transform  of  the  autocorrelation 

^^Kennedy,  Robert  S., Fading  Dispersive  Communication  Channels,  (New 
York:  Wiley-Interscience,  1969),  pp.  14. 

^^Ibid,  pp.l5. 
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Figure  4.29 


The  estimated  autocorrelation  function  shown  in  4.29  produces  the  frequency  response 
shown  in  figure  4.30  which  may  also  be  called  the  scattering  function,  s(co).  The  Doppler 
spread,  Bd,  of  the  channel  is  defined  in  Proakis'*^  as  the  range  of  values  over  which  the 
scattering  function  is  non  zero.  For  the  demonstrated  scattering  function  the  Doppler 


^^John  G.  Proakis,  Digital  Communications  (New  York:  McGraw-Hill,  1983), 
pp.  461. 
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spread  is  0.2  Hz.  The  coherence  time  from  Proakis^^,  (At)c“l/Bd,  is  5  seconds.  This 
frequency  spread  is  very  small  when  compared  to  the  usual  telemetry  transmission  signals. 

TTiis  response  looks  roughly  like  the  frequency  response  of  a  damped  oscillation. 
An  over  lay  of  this  response  is  shown  in  figure  4.31.  Hie  scattering  function  could 
roughly  be  modeled  with  a  damped  oscillation  response  with  a  mean  frequency  of  0.007 
Hz  and  a  3  dB  bandwidth  of  0.04  Hz.  The  transfer  function  for  this  model  is: 


s(a))  = 


_ 1 _ 

3.  33  -  i2.3(0  +  16.9(d" 


(35) 


This  appears  to  provide  a  useful  model  for  analysis.  The  agreement  is  within  the  variance 
for  most  frequencies  below  0. 13  Hz.  As  the  frequency  increases  the  noise  of  the  system 
begins  to  dominate  and  the  agreement  becomes  very  poor.  Figure  4.33  shows  the 
exponential  decay  is  present  in  higher  frequencies  for  the  model  and  the  estimated  scattering 
functions. 


44ibid. 
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Normalized  Response 


Scattering  Function  With  Damped  Oscillatitm 


Decibels 


6' 


Chapter  5 


Conclusions  and  Recommendations 

5.1  Conclusions  Concerning  Application  of  Results 

The  goal  of  this  thesis  has  been  to  look  at  the  effect  of  a  buoyant  plume  on  an 
acoustic  signal.  The  plume  was  looked  at  in  both  its  laminar  and,  at  least  partially, 
turbulent  forms.  Several  tools  were  used  to  perform  these  observations.  The  use  of  the 
Bom  and  Rytov  approximations  as  models  for  the  laminar  plume  was  investigated.  The 
scintillation  index  and  the  modulation  frequeiKy  spectrum  are  investigated  as  robust 
measures  of  turbulence. 

Specific  results  of  this  investigation  included:  the  verification  of  linear  spreading  in 
a  laminar  plume;  observation  of  a  x2  of  order  2  distribution  in  the  frequency  domain  of  the 
magnitude  square  response;  measurement  of  the  Doppler  spreading  due  to  the  plume  of  0.2 
Hz;  and  observation  of  a  damped  oscillation  frequency  response  model  for  the  scattering 
function. 

The  Bom  approximation  was  determined  to  be  a  good  measure  of  the  magnitude  of 
tlie  received  signal  but  did  not  follow  the  phase  changes  encountered.  The  Rytov 
approximation,  at  least  to  the  first  order,  is  not  as  good  as  the  Bom  approximation  in 
estimation  of  the  magnitude  response.  However  the  Rytov  approximation  provides  a  good 
estimation  of  the  phase  response.  Because  of  the  phase  response  and  the  small  absolute 
magnitude  of  the  error  in  magnitude,  the  Rytov  approximation  is  the  preferred  method  for 
modeling  this  system. 

The  frequency  distribution  of  the  modulation  of  the  received  amplitudes  indicate  the 
nature  of  the  plume.  As  the  plume  becomes  more  turbulent  the  spectmm  more  closely 
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approaches  a  Gaussian  distribution.  As  the  plume  develops  in  turbulence  the  spectrum 
initially  rises  at  the  higher  frequencies  and  at  the  fully  turbulent  state  the  lower  frequencies 
become  dominant.  This  seems  to  indicate  that  the  plume  at  first  creates  smaller  eddies,  and 
as  the  plume  becomes  more  fully  turbulent  the  larger  eddies  form  and  become  prevalent. 

From  the  experiments,  the  scintillation  index  proved  to  be  the  most  robust  estimate 
of  the  turbulence  in  the  plume.  It  provided  a  measure  of  the  transition  to  a  turbulent  plume 
in  both  the  well  controlled  fresh  water  experiment  and  the  large  salt  water  tank  experiment 
in  the  presence  of  uncontrolled  factors.  The  plume  should  be  interrogated  by  a  sufficiently 
strong  signal  and  the  pulse  train  should  be  of  sufficient  length  to  ensure  the  ability  to 
discern  the  signal  from  the  particulatly  scattered  signal.  In  addition  the  pulse  train  need  only 
be  of  a  length  to  allow  for  the  necessary  number  of  cycles  for  stable  statistical  properties. 
Hopefully,  the  length  of  the  pulse  train  will  be  significantly  less  than  the  period  of  the 
internal  waves  in  the  media  through  which  the  sampling  occurs. 

The  tools  used  in  this  thesis  provide  a  range  of  useful  methods.  Usir  ^  the  Rylov 
approximation  and  the  magnitude  and  phase  responses  of  the  signal  in  a  controlled 
experiment  can  provide  great  insight  into  the  plume  parameters.  The  frequency  response 
and  the  scintillation  index  provide  insight  into  the  controlled  experiment  as  well,  but  really 
offer  the  most  robust  tools  for  observing  a  plume  under  adverse  conditions. 


5.2  Recommendations  for  Further  Research 

The  overall  conclusion  of  the  research  conducted  is  that  there  needs  to  be  more 
research.  The  experiments  discussed  above  touched  on  many  aspects  of  the  thermal  plume 
without  fully  delving  into  any  single  area.  Several  experiments  immediately  suggest 
themselves. 
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A  detailed  study  of  the  plume  at  various  heights  using  multiple  sources  and 
receivers  simultaneously  on  the  same  plume.  TTiis  should  be  followed  by  a  detailed 
temperature  probe  investigation  of  the  plume  at  different  heights.  This  would  provide 
specific  constants,  similar  to  those  available  in  the  atmospheric  modeling  (eqns.  8  and  9), 
of  buoyant  plumes.  This  would  allow  specific  modeling  of  the  plume  based  on  the  acoustic 
interrogation.  With  limited  data  under  controlled  conditions  the  linear  spreading  of  the 
laminar  plume  was  verified.  More  turbulent  data  under  a  controlled  experiment  should 
verify  the  -5/3  spreading  of  the  turbulent  plume. 

An  experiment  placing  several  sources  and  receivers  around  a  buoyant  plume  at  sea 
with  a  set  geometric  orientation  to  provide  remote  sampling  of  the  plume  would  use  the 
observed  results  of  this  experiment.  By  using  a  pulsed  signal,  such  as  used  in  this 
experiment,  discrete  samples  of  the  plume  can  be  taken  and  the  robustness  of  the 
scintillation  index  as  a  measure  of  turbulence  of  a  plume  in  a  natural  environment  can  be 
assessed. 

Tlie  submarine  spring  plume  is  a  common  buoyant  plume  available  for  monitoring. 
Colboume  and  Hay'^^  conducted  a  study  of  one  such  plume  using  a  submersible  and 
acoustic  backscattering.  The  source  was  roughly  5  meters  in  diameter  with  vertical 
velocities  near  30  cm/s.  Salinity  measured  at  one  meter  varied  up  to  5  ppt  while  temperature 
varied  only  about  0.2°  C.  This  leads  to  a  peak  velocity  differential  in  the  plume  of  5.2  m/s 
due  to  the  dominant  effect  of  the  change  in  salinity  for  these  values.  A  change  of  5.2  m/s 
over  5  meters  produces  a  very  significant  effect  on  the  propagated  signal  and  should  be 
easily  measurable  (the  experiment  in  this  paper  produced  roughly  15  m/s  change  over  2 
cm). 


^^Eugene  B.  Colbourne  and  Alex  E.  Hay,  "An  Acoustic  Remote  Sensing  and 
Submersible  study  of  an  Arctic  Submarine  Spring  Plume, "Journal  of 
Geophysical  Research,  95,  No.  C8  {15  August  1990),  13219-13234. 
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In  recent  years  another  source,  the  hydrothermal  vent  fields,  has  been  subject  to 
considerable  research.  Little,  Stolzenbach,  and  Von  Herzen'*^  used  DSRA/^  Alvin  to  obtain 
measurements  in  a  vent  field  and  obtained  vertical  velocities  near  5-lOcm/s  and  temperature 
peaks  slightly  greater  than  2°C  at  24  meters  above  the  field.  Similarly,  Speer  and  Rona”*^ 
arrived  at  similar  values  in  their  models  of  hydrothermal  plumes.  These  fields  are  highly 
turbulent  and  extend  over  100x100  meter  areas.  On  a  smaller  scale,  single  hydrothermal 
openings  are  on  the  order  of  2-8  cm.  Temperature  variations  closer  to  the  vents  are  from 
30-420°  C  for  hot  smokers  and  5-30°  C  for  diffuse  flow  from  cracks.'**  There  currently 
exists  no  remote  sensing  systems  monitoring  any  of  these  fields  on  a  continuous  basis. 

The  possibility  of  using  forward  scattering  across  one  of  these  vent  fields  as  an 
economical  means  of  monitoring  the  activity  of  the  field  offers  an  obvious  application  of 
this  research.  A  small  network  of  several  transducers  and  hydrophones  would  offer 
different  paths  for  monitoring  a  large  area. 

Finally,  the  nature  of  buoyant  plumes  due  to  variation  in  the  layered  media  by  a 
horizontal  body  motion  can  be  investigated  using  this  line  of  sight  method.  The  time  series 
plots  of  signal  response  obtained  from  an  experiment  involving  the  insonification  of  the 
wake  left  by  22.9  cm  diameter  by  122  cm  slender  body^®  appears  to  be  similar  to  the  type 
of  responses  obtained  in  this  experiment.  Analysis  by  both  the  scintillation  index  and  the 
frequency  transformation  of  the  time  series  discussed  above  may  reveal  significant 
information.  This  type  of  frequency  analysis  has  not  been  done  in  the  available  literature. 
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